Static Fundamental Solutions of Einstein Equations 
and Superposition Principle in Relativistic Gravity 



P. P. Fiziev* 



Abstract 



We show that Einstein equations are compatible with the presence of massive 
point particle idealization and find the corresponding two parameter family of so- 
lutions. They are complete defined by the bare mechanical mass M > and the 
Keplerian mass m > (m < M) of the point source of gravity. The global analytical 
properties of these solutions in the complex plane define a unique preferable radial 
variable of the one particle problem. 

These new solutions are fundamental solutions of the quasi-linear Einstein equa- 
tions. We introduce and discuss a novel nonlinear superposition principle for so- 
lutions of Einstein equations and discover the basic role of the relativistic analog 
of the Newton gravitational potential. For the relativistic potential we introduce 
a simple quasi-linear superposition principle as a new physical requirement for the 
initial conditions for Einstein equations, thus justifying the instant gravistatic case 
for N particle system. 

This superposition principle allows us to sketch a new theory of the gravitational 
mass defect. In it a specific Mach-like principle for the Keplerian mass m is valid, 
i.e. it depends on the mass distribution in the universe, in contrast to the bare mass 
M, which remains a true constant. Several basic examples both of discrete and of 
continuous mass distributions are considered. 

1 Introduction 

1.1 Static Fundamental Solution and Superposition Principle in 
Newton Theory of Gravity 

The notion of a static fundamental solution of a classical field equation appeared at first 
in Newton theory of gravity [I]. Such a solution solves the Poisson equation with source 
term, proportional to Dirac 3D 5-function: 



Here and further on for simplification of the formulas we are using units in which the 
Newton gravitational constant Q Newton — i an d velocity of light c = 1. 
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The static fundamental solution does not depend on arbitrary functions, or additional 
constants. It is unambiguously fixed, among all solutions of Eq. (II. ip . if we require this 
solution to tend to zero at infinite distances. This unique solution describes the potential 
of the newtonian gravitational field 
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r 

which is created by a classical point particle of gravitational (Keplerian) mass m, placed 
in Euclidean 3D space at the origin of coordinate system r = 0. Here r = |r| > 0. 

Analogous solutions are well known in the problem of static point source of electric field 
in Maxwell electrodynamics, as well. A proper generalization of the notion of fundamental 
solution for hyperbolic partial differential equations can be find, for example, in [2] and 
in the references therein. 

According to well known mathematical results, the solutions ( II. 2p describe, too, the 
static field in vacuum, outside sources of finite dimension, assuming spherical symmetry 
of the corresponding distribution fi(r) = fi(r) of gravitational mass, or electric charge. 

The fundamental role of the solutions (11.21) in Newton gravistatics is substantiate by 
the superposition principle, according to which in these linear theories the field of any 
aggregate of matter can be obtained as a sum of the fields of its constituent matter points 
at positions (or r'). For example, in Newton gravistatics 
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in the case of a set of N discrete massive points, and 
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in the case of a continuous distribution of Keplerian mass fi Kepler (r). 

In a more general form the superposition principle in the newtonian gravity can be 
expressed as 

Proposition 1: If /i/(r) and /i/j(r) are two mass distributions, which create gravita- 
tional fields with corresponding potentials ipf emton (r) and (pfj wton (r), then the potential of 
the field, created by mass distribution yu(r) = /i/(r) + fin(r) is 
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<P ( r ) = ^/ W + V// ( r )- (1-5) 



1.2 Massive Point Particle in General Relativity 
1.2.1 The Schwarzschild Solution 

An attempt to solve the point-mass problem in general relativity (GR) was made at 
first ninety years ago, as early as in the pioneering article by Schwarzschild [3] and its 
subsequent modifications [HE]. 

Today it is well known that inconsistences arise when we look at Schwarzschild solution 
as the space-time arising from localized point mass singularity [B]. Actually, the well 
known Schwarzschild metric in Hilbert gauge: 



solves the vacuum Einstein equations G^ = in the spherically symmetric static case. It 
possesses an event horizon at p = pa = 2m and a strong hidden singularity at p = 0. 
This solution describes a completely empty space-time with removed point p = and 
nontrivial non- Euclidean topology. Indeed, in the Weyl's isotropic coordinates with radial 

variable r w = | — m + a/ p{p — 2m) j one can easily see that this solution describes a 

two sheeted space-time with a flat asymptotic at rw = and rw = oo connected by a 
specific bridge. Then the nonzero Keplerian mass m appears in the solution due to the 
nontrivial topology of the space-time, i.e. in the spirit of Einstein-Rosen-Misner- Wheeler 
geometrodynamics |7J, as described by the sentence "mass, without mass". 

The singularity at p = is not related with a massive point particle with proper bare 
mass M and mechanical action Am = —M J ds. Indeed, as we shall see in more detail 
below, the solution (jl.fip does not solve the Einstein equations (EE) 

= 8tcT? (1.7) 

in presence of matter with stress-energy tensor T>? ~ M^(r). Here M5 (3) (r) describes 
the mass distribution of the point particle with proper bare mass M. 

In his pioneering article Schwarzschild has used another radial variable r, which defers 
essentially from variable p, i.e., his choice a radial gauge for the spherically symmetric 
static metric 

ds 2 = g tt (r ) 2 dt 2 + g rr (r ) 2 dr 2 - p(r ) 2 (d6 2 + sin 2 9 d(f) 2 ) (1.8) 

is different from the Hilbert's one (11.61) . 

Borrowing from the Minkowskian flat space-time the gauge condition 

\ 4 g\:=det\\g IB ,{t,T)\\ = l 

(which takes place there in Cartesian coordinates), Schwarzschild was able to fix the three 
unknown functions in the form: 

I 2m 

p(r) = {jr* + p% > 0, g tt {r) = 1 - — > 0, and g rr (r) = -p 1 \r) 2 / g tt (r) < 0. 

This solution of EE has no event horizon. Its peculiar feature is that it describes a 
point like object of Keplerian mass m > 0, zero radius, zero volume, but nonzero area 
A PG = 4tt p 2 G > 0. These unusual properties of the original Schwarzschild solution have 
been discussed by Brillouin [8] as early as in 1923. 

At present the original Schwarzschild geometry and other similar geometries of space- 
time are widely ignored in GR. A main stream of articles in the last 40 years is strongly 
limited to consideration of the black hole interpretation of the Hilbert form (II. 6p of vac- 
uum Schwarzschild solution and its generalizations. In addition, essential features of the 
original Schwarzschild solution are not reproduced in the most of modern literature on 
this subject and remain hardly known. 



1.2.2 On the Choice of Radial Variable in the GR Massive Point Particle 
Problem 

One of the basic problems in the description of single massive point particle as a source 
of gravitational field in GR is the choice of proper radial variable r. 
The quantity p > has a clear geometrical and physical meaning: 

i) It is well known that p defines the area A p = 4np 2 of a centered at r = sphere 
with "area radius" p and the length of the big circle on it l p = 2np. Thus we see that the 
quantity p has a well defined geometrical meaning and is a gauge invariant notion. 

ii) The coordinate p measures the curvature ~ 1/p 2 of the 2D-manifolds (2D-spheres) 
in 3D Riemannian space, which are invariant under rotations around the center of spherical 
symmetry It measures, too, the curvature of the 4D pseudo-Riemannian space-time: 
4 R = 4 R{p) and of the corresponding 3D-space: 3 R = 3 R{p) in the spherically symmetric 
case, inside the matter source of finite dimension. Hence the name "curvature radius". 

iii) From physical point of view one can refer to p as an optical "luminosity distance", 
because the luminosity L of distant physical objects is reciprocal to A p : L ~ 1/ p 2 . 

In contrast, the physical and geometrical meaning of the coordinate r is not defined 
by the spherical symmetry of the problem and is unknown a priori [9] . Its choice has been 
discussed from physical point of view by Eddington as early as in [10J. His conclusion 
was that all admissible variables r are practically equivalent at distances r >> pc, since 
under suitable coherent choice of their scales we have p/r — > 1 when r — > oo. 

The following assumptions about the mathematical properties of the radial variable r 
of the single point particle problem seem to be natural from physical point of view: 

i) Its value r = is to correspond to the center of the symmetry, where one must place 
the physical source of the gravitational field - the massive point particle. 

ii) The radial variable is to vary in the semi-bounded physical interval r e [0, oo). 

iii) The luminosity variable p(r) is to increase monotonically to infinity in this inter- 
val, together with radial ones, i.e., dp/dr > 0, and, in addition, one has to impose the 
Eddington condition: 

lim M = 1 

r—>QO r 

iv) The infinite value r = oo of the radial variable is to be prescribed to the boundary 
of the asymptotically flat domain of space-time. 

v) There must not exist non-physical singularities of the solution of EE in the whole 
compactified complex domain C of the radial variable r G C. 

In the present article we show that these physical requirements define in a unique way 
the radial variable r of the problem at hand, thus solving the corresponding uniformization 
problem in the case of point particle source of gravity in GR. 

1.2.3 Some Remarks on the Massive Point Particle Idealization in GR 

A clear physical motivation for consideration of massive point particle sources of gravi- 
tational field in GR, both electrically neutral and charged ones, can be found in 1962-63 
Feynman lectures on gravity [llj . The energy momentum tensor of a point particle has 
been used in the excellent textbooks by Landau & Lifschiz and by Weynberg [1] as a tool 
for treatment of many particle systems in GR. In spite of this fact the single particle case 
is still an open problem in GR. 



(1.9) 



Moreover, at present the vast majority of relativists do not accept the consideration 
of point particles in GR, assuming that it is an idealization, which is incompatible with 
EE [12]. There are different reasons: 

i) Some doubts about consistence of the theory of mathematical distributions (like 3D 
Dirac 5-function 5®(r)) [13] with the obviously nonlinear character of EE [T3] . 

The formal mathematical problems, which emerge when one attempts to work with 
distributions in EE were successfully advanced in the last decade using Colombeau's 
theory of generalized functions [15]. Unfortunately, the published results on the point 
particle problem in GR, based on this approach, are physically incorrect (see Section 
2.3.5). 

ii) The clear understanding that an infinite concentration of energy in a single space- 
point will change drastically the geometry of the GR-Riemannian space- time M^ 1 ' 3 )^} 
in a small vicinity of the world-line of this point; 

iii) Some attempts to neglect the role of classical description of matter in GR, replacing 
it by classical field description, or by quantum field description, according to the so called 
"third approach" by Einstein, Wheeler and many others, see [7J and references therein; 

iv) The absence of understanding of necessity to use fundamental solutions of EE. 
These were unknown up to recently in GR, but may turn to be useful mathematical tool. 

v) The absence of a general non-linear superposition principle for EE, which is to 
correspond to the linear superposition principle (11.31) - (11. 51) in Newton theory of gravity. 
Note that a specific kind of superposition principle for initial conditions of black hole 
solutions is well known [TJ, but the general problem and other specific cases are still not 
studied, to the best of our knowledge. 

On the other hand, it is obvious that in Nature the very distant stars look like "points" 
of finite mass and finite luminosity. This fact has a proper mathematical description in 
the language of mathematical distributions in the Newton theory of gravity, but still not 
in GR. 

In spite of the absence of proper description of the massive point particles, in the 
practical relativistic celestial mechanics, for example, in the calculations of the solar- 
system trajectories of the space crafts, even the Sun is considered often as a massive 
point source of gravity. 

A formal mathematical problem is to find the corresponding correct treatment of such 
objects in GR, but up to recently no reasonable approach was known. Unfortunately, the 
most of the existing formal attempts to solve the point particle problem in GR do not 
take into account one essential physical difference between the GR and the Newtonian 
description of the massive objects. It is well known [TJ, that any body in GR has two 
different masses: the Kepelrian one m, as seen from distant observer, and the proper 
(bare) mass M > m, which is the sum of the masses of its constituents, when placed 
at infinite distances between them, i.e., with gravitational interaction - turned off. The 
difference M — m, or the ratio g = m/M describe the gravitational defect of mass. One 
must include properly this specific feature of the relativistic theory of gravity in the GR- 
point-particle model. To the best of our knowledge, such attempts ware not made up to 
recently, with the only exception - [16J. 

In the present article we show that a correct mathematical solutions of EE with 5^ (r) 
term in the rhs do exist. Such solutions describe a two parameter family of analytical 
space-times M^ 1 ' 3 ^^^} with a specific strong singularity at the place of the massive point 
source with bare mechanical mass M > and Keplerian mass m < M. 



The price, one has to pay for this enlargement of the standard GR framework, is: 

i) To accept the unusual geometry of the space-time around the matter point with 
infinite concentration of energy in it. 

As we have stressed already, this geometry was introduced in GR for the first time 
actually in the original Schwarzschild article [3] and has been discussed by Brillouin [8]. 
The unusual geometry is essentially different from the geometry around the space-time 
points with finite energy density in them. 

The global properties of the space-time manifolds, generated by massive point source, 
are essentially different from the ones of the most popular, at present, Hilbert form (11.61) of 
the original Schwarzschild solution, which describes an empty space-time with nontrivial 
topology in the spirit of the Einstein- Rosen- Wheeler geometrodynamics [?]. 

ii) To allow consideration of metrics, whose coefficients are not a C 3 - smooth functions. 
Indeed, to reproduce the 8^(r) term in the rhs of EE, the metric tensor, and/or its 
derivatives, related to the geometry of the Riemannian space-time, must have a definite 
singularities (discontinuities) at the place of the point source of gravity [16] . 

iii) To replace the mathematical theory of the real smooth manifolds, which is in 
current use in GR, with the theory of the analytical manifolds with proper singular points, 
considering the whole complex domain of the space-time variables. 

1.2.4 The Nonlinear Superposition Principle in GR 

In addition, here we formulate a GK-nonlinear superposition principle, analogous to the 
Newton one, described in Proposition 1. 

Let's consider for simplicity only the case of asymptotically flat space times, which cor- 
respond to energy- momentum stress tensors Tjf{x) with compact support in M^'^jg^}, 
i.e. let's focus our attention on so called "island universes". It is well known that after a 
proper fixing of the gauge, the boundary conditions at infinity define the solutions of EE 
(11.71) in a unique way, see for example [TJ1 [1] and the references therein. 

Proposition 2: Let g^ v (x)j and g^^jj are two metrics, which correspond via EE 
( | 1. 7| ) to two energy-momentum tensor distributions Tjf(x)j and T^(x) II of compact sup- 
ports. Then the metric g^ix) of the GR gravitational field, created by energy-momentum 
tensor distribution T£{x) = T^(x)j + Tff{x) n via EE ( [i. 7[ ) is uniquely defined by the two 
metrics 5 f A1 i/(x) / and g /JiV (x) II . 

Thus we obtain the unambiguous correspondence 

{9^(x) I ,g flu (x) II } i-> g^(x) (1.10) 

and the metric g^ y {x) deserves to be called a nonlinear superposition of the metrics g fll/ (x) I 
and g^{x) ir 

The essence of the proof of the existence of such nonlinear superposition principle 
in GR is in the simple note that the support of the distribution Tjf(x) = Tjf{x) I + 
Tjf(x)jj will be certainly compact, if both Tj^{x) I and Tjf(x) H have compact supports 
in Ma 1 ' 3 )-^^}. Then the EE (11.71) with Tjf{x) in rhs, supplied with asymptotically flat 
space-time boundary conditions, will have an unique solution g^ u {x), which corresponds 
to the metrics g flv (x) I and g flv (x) II and ought to be named their (nonlinear) superposition 
due to obvious physical reasons. It is clear that g^u(x) is a very complicated functional of 
the metrics g^(x)j and g fJiU (x) II . The problem of reconstruction of g^ u (x) in (11.101) . using 



two given metrics g fJiV (x) I and g, J ,v{x) II is highly nontrivial. We shall use the symbol (S) 
to denote the composition (11.101) of the two metrics in the form 



9nu{x) = g, a ,{x) I ®g IJa ,(x) II . 



(1.11) 



It is clear that by construction this new operation on the metrics is symmetric and asso- 
ciative: 



These properties are an immediate consequences of the corresponding properties of the 
summation of energy- momentum tensors of compact supports in M^' 3 \ assuming that 
we are considering space-times with a fixed flat-geometry-boundary-conditions at space 
infinity. 

Finally, we can define correctly the superposition of an arbitrary number N of met- 
rics, obeying the same boundary conditions at infinity, i.e. we can introduce a multiple 
superposition operation: 



To some extend the novel principle (11.111) . HI . 13|) is unexpected, and certainly much 
more complicated than the simple linear superposition principle (II. 5p in Newton gravity. 

Indeed, in GR we have a very specific physical situation. It is clear that even if the 
metrics g IJ , l ,(x) I and g lxv {x) II are static, in general case their superposition g^(x) is not 
a static metric. It contains the whole GR dynamics, including the possible radiation of 
gravitational waves, due to the gravitational interaction between the physical sources of 
the metrics g^u^j and g llu (x) II . In contrast to the situation in electrodynamics, where we 
can introduce non-electrodynamical forces and stresses to keep the composite aggregate of 
charges in a static state without introducing new terms in the Maxwell equations, in GR 
any additional interactions, introduced for the same purpose, will have a nonzero energy- 
momentum tensor T^(x) additional , which enters the rhs of EE and changes the space-time 
geometry and the very problem. As a result we see that without introducing of non- 
gravitational interactions between particles in GR we have only a unique (whole-time) 
static case - the single point particle problem (see, for example, Fock in pQ). 

Nevertheless, as we shall show in the present article, considering just the instant static 
case, one can introduce a simple quasi-linear superposition principle for static fundamental 
solutions in GR. It reveals the role of static fundamental solutions of EE in GR, which 
is much like the role of corresponding fundamental solutions in linear field theories like 
Newton gravity and Maxwell electrostatics, if one considers only a single 3D space-time 
surface t = [T7] . 

The resolution of the relativistic gravistatic problem requires the solution only of the 
well known suitable form of it -EE. It does not contain second derivatives of metric with 
respect to the time variable t. Let us consider a space-time, which is a solution of EE. 
The 3D curvature of arbitrary 3D space-like surface in it obeys the well known basic 
equation 




(1.12) 



9^u{x) = ® A=I g^ v {x) A := g (tt/ {x) I ® ... ®g^(x) 




{3) R + K 2 = l&Kfi. 



(1.14) 



Here Ki = K 2 — K^K 1 ^ (where K = g^Kij) is the exterior curvature of the 3D surface 
and \x is the relativistic density of mass distribution. 

If the last equation is fulfilled at some time instant t, as a consequence of the EE it 
will be fulfilled at all time instants t G (—00, 00), for which the problem is well defined. 
There exists an inverse theorem, too: The whole system of EE, which governs dynamics 
in GR, may be derived by the requirement to have the above relation (11.141) co-variantly 
valid at all time instants |18j . 

According to articles [T7j, one can define the relativistic gravistatics as a description 
of time-symmetric initial value problem for EE in proper coordinates, as well. To see 
this, it is enough to know that choosing appropriate coordinates outside the 3D surface 
t = one obtains for the coefficients of the second fundamental form of this surface 
Kij = —^dtgij. Then one defines the instantaneous-static solutions of EE, possessing a 
3D space-like surface t = on which = 0. The coordinate independent way to this 
definition implies existence of isometry of space-time: t — > —t, r — > r [17J. Then the basic 
equation reduces to the following simple form, valid at time instant t = 0: 

(3 lR = 16tt/x. (1.15) 

This equation is not a dynamical equation, but just a constraint on the initial con- 
ditions - a specific relativistic constraint equation (RCE). As a result of relativistic dy- 
namics, governed by the other EE, the RCE will be automatically fulfilled for any time 
t, if it will be valid at time instant t = [T71 [T]. Hence, the time t is a simple auxiliary 
parameter in the RCE and in its solutions. Thus we see that it is enough to solve the 
RCE only at the initial time instant t = 0, i.e., it is enough to solve the Eq. (jl,15p . 

From pure mathematical point of view any of the solutions of Eq. (11.151) may be 
considered as an initial condition of a proper initial value problem for EE. One of the 
basic purposes of present article is to find the physically meaningful solutions of RCE 
among the whole variety of its possible mathematical solutions. These physical initial 
conditions were not known until now. Their discovery calls for reconsideration of many 
well studied problems in GR, including the gravitational collapse problem. 

We present here an instant static solutions of the RCE with singularities, which cor- 
respond to presence of arbitrary number of massive matter points, both of discrete or 
of continuous distribution. These solutions define the physically meaningful initial con- 
ditions for EE, which describe the real matter, made of massive point particles. They 
present a very special class in the variety of all initial condition, which are admissible 
from pure mathematical point of view. 

To obtain the non-stationary gravitational field of moving matter point sources, the 
instant solutions of RCE can be modified in a manner, which is well known from relativistic 
electrodynamics. Of course, the whole problem is a highly complicated and still not 
solved. One can hope that in GR a procedure, which is analogous to the introduction of 
Lieanard-Wiechert potentials may take place. Here we shall stress that in electrodynamics 
for this purpose the static Colomb potential is in use. We believe that in a similar way 
our fundamental solutions may turn to be the key tool for the correct treatment of the 
nonlinear GR dynamics of many-particle systems, which is still an open problem. 

As already mentioned, in GR, due to gravitational mass defect, we have to distinguish 
two different masses of every body - the Keplerian gravitational mass m and the proper 
bare mass M > m of the body pQ. Under some additional natural assumptions our 
new quasi-linear superposition principle for the fundamental solution of EE yields a novel 



theory of the relativistic gravitational mass defect of systems of discrete matter points 
and composite bodies of continuous mass distribution. It is based on a specific integral 
equation for the relativistic gravitational potential, derived for the first time here. It turns 
out that the GR mass defect is governed completely by the RCE, with time t, playing 
the role of an auxiliary parameter in its solutions. Thus, our quasi-linear superposition 
principle for the fundamental solution of EE has a basic impact on the relativistic theory 
of gravitational defect of mass. We give here for the first time a number of solutions of 
the integral equation for the relativistic mass defect. These describe some basic physical 
problems: two-point-particle problem, some special cases of three and four point particle 
problems, the general properties of the iV-point-particle problem, as well as some basic 
examples of continuous mass distribution: homogeneous massive circles, spheres and balls. 

Further important physical consequences, which can be derived using the new funda- 
mental solutions of EE and the corresponding nonlinear superposition principle for them, 
will be considered elsewhere. 

2 The Mathematical Problem of Single Point Parti- 
cle in GR 

2.1 The Total Action and Introduction of Coordinates 

Let us suppose that in the whole universe there exist only a single massive point particle 
of bare mechanical mass M, and that it creates its own gravitational field according to 
the laws of GR. This problem is described by the total action A tot =Aqr -\-Am- The first 
term describes the action of the gravitational field, created by the single particle. The 
second term adds to the total action the pure matter (mechanical) action of the massive 
particle. Thus in GR the total action acquires the well known explicit form: 



In the rest frame of the point particle both the total action (12 .ip and the formed 
by this particle GR space-time manifold M^ 1 ' 3 ^^^} have an obvious group of symmetry 
SO (3) x T ( (l). As a result, the problem can be reduced not only on the orbits of the 
group 5*0(3), i.e. on the 2D quotient space 

M (i,i) = M^/SO(3), with natural global 
coordinates t and r, but even on the orbits of the whole group S0(3) x T t (l), i.e. on the 
ID quotient space M« = M^/(S0(3) x T t (l)), with some natural radial coordinate r. 

To be able to use some coordinates x = in the Riemannian space- time M^ig^} 
of the point particle problem, one actually presupposes to have a flat Minkowskian space- 
time E^ 1 ' 3 )-^^}, endowed with the same coordinates. For example, one assumes to bor- 
row the Cartesian coordinates: {£, r}, or the spherical ones: {t, r, 8,(p} from the flat 
Minkowskian space-time for the use in the Riemannian space-time. Thus we have at our 
disposal simultaneously a flat metric 77^ (x), and a Riemannian metric g^ u (x), expressed 
in the same coordinates. 

One of the basic results of present article is that the auxiliary flat space- time E^ 3 )^} 
plays much more profound role in the problem at hand, than the usually expected formal 
one. In particular, it turns out that the real geometrical points of the two space-times: 
E^^r^l/Wo, with the world line W of one point (the origin r = 0) removed, and 




(2.1) 



^ ~~ ^ e space-time of the GR massive point particle problem, are in one-to-one 

correspondence. In particular, the last has the same topology as the first one. Note 
that the Galilean space-time of a single point particle problem in Newtonian gravity has 
precisely the same topology as the one of the manifold E^'^jr^j/Wo. Thus, in the real 
domain of variables, the space-time in the Newton gravity and in the GR, formed by a 
single massive point particle, have the same topology. This observation makes it clear that 
even because of pure topological reasons the black hole solutions of EE are not compatible 
with the matter point sources of gravity, sice they have a different topology. 

In its proper frame the single massive point particle, placed at the origin of the stan- 
dard spherical coordinate system in the 3D Riemannian space M^lg^-} C M^'^jp^} 
yields the familiar static metric (11. 8ft with three unknown functions gu{f) > 0, g r r{r) < 0, 
and p(r) > of the radial variable r > [TJ. The variable r is not defined by the 50(3) 
symmetry of the problem, nor by its global-time translation invariance with respect to 
the group T t (l). From geometrical point of view the choice of the function p(r) fixes the 
imbedding of the quotient space 

M (i) = m( 3 V50(3) into the 3D space M®. We assume 
that by definition the value r = of the radial variable r corresponds to the center of 
spherical symmetry, C. There the massive matter point is placed. We also accept other 
assumptions about the mathematical properties of the radial variable r, listed in the 
Subsection 1.2.2 of the Introduction. 



2.2 On the Role of the Gauge Fixing in the Massive Point Prob- 
lem 

General relativity is a gauge theory. The fixing of the gauge in GR is described by a 
proper choice of the quantities 



I\ = )=g^d x (^\g~\g x » 



in the 4D d'Alembert operator g^V ^ V ' v = g^ v (d^d v — Y^dS) [TJ, and actually is a fixing 
of the coordinates. In our problem the choice of spherical coordinates and static metric 
dictates the form of three of the quantities f^: f t = 0, f g = —cot 9, = 0, but the 
function p(r) and, equivalently, the form of the quantity 



r r = (In 



gitp 2 



are still not fixed. Here and further on, the prime denotes differentiation with respect 
to the variable r. We refer to the freedom of choice of the function p(r) as a rho-gauge 
freedom in a broad sense, and to the choice of the p(r) function as a rho-gauge fixing. 

At first glance the function p(r) may be chosen in quite arbitrary way, thus fixing the 
remaining (radial) gauge freedom of the problem - the only one, which is not fixed by 
symmetry reasons. We show that choosing a definite class of functions p(r) one can solve 
correctly the EE ( 11. 71) with stress-energy tensor 

A(3) f r ^| 

2? = M<J< 3 VK^ = M-==U=<^<$° (2.2) 



which describes a massive point source with bare mass M at rest in stationary and static 
coordinates. It may seem strange that for solving this problem, one needs to fix the class 
of coordinates by a proper choice of the radial gauge. As we shall see, the choice of the 
admissible class of radial coordinates r is a consequence of the boundary conditions. In the 
problem at hand these conditions are masked in 3D Dirac (^-function in (12.21) . It describes 
in a formal mathematical way the properties of source of gravity and its boundary. 
The following comments throw an additional light on this delicate issue: 

1. The strong believe in the independence of the GR results on the choice of coordi- 
nates x in the space-time M^ 1,3 ) {g^ v (x)} predisposes us to a somewhat light-head attitude 
towards the choice of the coordinates for a given specific problem. Indeed, it is obvious 
that physical results of any theory must not depend on the choice of the variables and, 
in particular, these results must be invariant under any admissible changes of the coor- 
dinates. This requirement is a basic principle in GR. It is fulfilled in any already fixed 
mathematical problem. 

2. Nevertheless, the change of the interpretation of the variables may change the 
formulation of the very mathematical problem and thus, the physical results. This can 
happen, because we are using the variables according to their meaning. For example, if 
we are considering the luminosity distance p as a radial variable of the problem, it seems 
natural to put the point source at the point p = 0. In general, we may obtain a physically 
different problem, if we are considering another variable r as a radial one. In this case we 
shall place the source at a different geometrical point r = 0, which now seems to be the 
natural position for the center C, but does not coincide with the previous one - p = 0. 
Imposing the same physical requirements, i.e. the same boundary conditions at different 
places in the space, we obviously will obtain different physical problems and results. Of 
course, as in any gauge theory, in GR there exist a classes of physically equivalent gauges. 
All gauges (coordinates) in such class yield the same physical results. The real problem 
is how to find the correct class of the gauges, proper for the given physical configuration. 

3. The relation between the two geometrical "points": p = and r = 0, and between 
the corresponding physical models of a point particle, strongly depends on the choice 
of the class of functions p(r), i.e. on the class of the radial gauges. Thus, applying 
the same physical requirements in essentially different "natural" variables, we arrive at 
different physical models, because we are solving EE under different boundary conditions, 
coded in corresponding 3D Dirac 5-functions in (12. 2p . One has to find a theoretical or an 
experimental reasons to resolve this essential ambiguity. 

4. The choice of the radial coordinate in the single point particle problem in GR needs 
a careful analysis. It is essential for the description of the very source of gravitational 
field, not for the description of the field in surrounding this source vacuum domain. A 
well known mathematical fact is that in the vicinity of a definite singular point of a 
mathematical functions one must use a definite special type of coordinates for adequate 
description of the character of the singularity, i.e., one is to solve the corresponding 
uniformization problem. 

5. The solutions of EE in essentially different coordinates have different singularities 
somewhere in the whole complex domain of the corresponding variables. The essentially 
different coordinates may be equivalent only locally - in the spirit of the widely used 
theory of smooth manifolds. One ought to make a reservation, speaking about "essentially 
different coordinates" , because there exist a coordinate changes, which alter only the place 
of the singularities of the solutions of EE in the complex domain of the variables, without 



varying the character and the number of these singularities. Such changes are precisely the 
linear ones and the fractional-linear ones. All other, more general coordinate changes, do 
not possess such property and yield essentially different coordinates in the whole complex 
domain. 

6. In our particular problem, according to Birkhoff theorem, the spherically symmetric 
solution with given Keplerian mass m is unique in the vacuum domain. The coordinates, 
which are essentially different somewhere else, may be locally equivalent in the vacuum 
domain. As a result, all local GR effects, like gravitational redshift, perihelion shift, 
deflection of light rays, time-delay of signals, etc., will have their standard exact values 
in static spherically symmetric gravitational field with given Keplerian mass m. These 
physical values do not depend on the admissible coordinate form of the solution. 

We will use this local gauge freedom in description of the gravitational field outside 
the source to reach an adequate mathematical modelling of the very point source. 

2.3 The Gravitational Field Equations and Their Solution 

2.3.1 The Vacuum Solution in an Arbitrary Radial Gauge 

The EE (ll.7p for our problem with metric (If .81) can be easily derived from the following 
form of the nonzero components of Einstein tensor: 

2 (£\- 3 (tf+a££E)+L (2.3a) 
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using, in addition, the corresponding components of the energy-momentum tensor (I2.2p 
of a single matter point. 

In particular, solving the EE for the static, spherically symmetric case in vacuum, one 
easily obtains the following most general solution: 

g tt (r) = 1 - ^- > 0, g rr (r) = - (p'(r)) 2 / g tt {r) < 0, 
p(r) 

p(r) - an arbitrary C 1 function. (2-4) 

It was derived for the first time already in the articles |19j . 

If one uses the Hilbert gauge Ph{ t ) = ^ in the EE with 8(r) term in the rhs, one easily 
reaches a contradiction [HI US]. Hence, now the question is as to how to choose the radial- 
gauge-fixing-f unction p(r), to be able to comply with the specific boundary conditions at 
r = 0, coded in the 5(r) term in the rhs of EE. In other words, we have to find a radial 
gauge, which makes the boundary problem for EE consistent with the presence of matter 
point source of gravity. 



2.3.2 Normal Coordinates for Static Spherically Symmetric Gravitational 
Field 



The expressions (I2.3P demonstrate a very important feature of EE: In spite of their non- 
linearity, which may yield doubts in the applicability of the theory of mathematical dis- 
tribution, EE are quasi-linear differential equations. After all, the higher (second) order 
derivatives of the unknown functions enter these equations linearly. This makes possible 
the usage of mathematical distributions [13] in the GR massive point particle problem 
in some specific coordinates and the usage of the Colombeau's theory of the generalized 
functions [15], hoppefully in all admissible coordinates. 

The fundamental solutions of EE were found for the first time in the articles [T6] . 
introducing (in a slightly different notations) proper normal field variables <p(r), <f2{r) 
and <f(r) according to the formulas 

g tt = exp(2<p), p = pexp(-<p + y? 2 ), g rr = - exp(-2<p + 4y? 2 - 2<p). (2.5) 

Here p = const > defines the scale of the luminosity variable. 

In the present article we develop a more general approach to the fundamental solutions 
of EE, which allows consideration of arbitrary number of massive point sources of gravity. 
We shall see that an essential ingredient of this approach remains the Fock conformal 
transformation of the space MS 3 \ see Fock in [lj. It arises naturally when one puts the 
gravitational action of the static spherically symmetric problem to a canonical form [16] . 
Indeed, after reduction of the Hilbert-Einstein action Aqr on the orbits of the group 
5*0(3) x T t (l), we arrive at one dimensional variational problem with "Lagrangian" 

t l/ WUgVO^ +vS ^-\ (2.6) 

The corresponding Euler-Lagrange equations for pure gravitational field in vacuum read: 
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V~9rr = 0, (2.7a) 
0, (2.7b) 



9ttV-9rr = 0, (2.7c) 



\J~9rr 

where the symbol " = " denotes a weak equality in the sense of theory of constrained 
dynamical systems. As a result of the rho-gauge freedom the field variable \J—g rr is not 
a true dynamical variable but rather plays the role of a (specific nonlinear) Lagrange 
multiplier, which is needed in a description of constrained dynamics. Its derivative with 
respect to the radial variable r does not enter the Lagrangian (12.61) . An advantage of such 
derivation of field equation is that it makes transparent this fact. Of course, the equations 
(12.71) are completely equivalent to the vacuum EE, considered in the previous subsection. 
As a result equations (12.71) are solved by the functions (12. 4p . 

Let us consider the formal 2D space M^ 2 ) of the field variables J~g~tt and p, endowed 



with the quadratic metric form . 2p dpdi^/gTt) H — dp 2 . It is easy to check that 

\/-9rr{p) V y'-arr(p) 



its Riemannian curvature tensor is zero. Hence, in this space one can introduce a normal 
field variables, transforming its 2D metric into canonical form. The above change of 
variables yields the corresponding diagonal form of the Lagrangian: 

L=\(f (-( P V) 2 + (pV 2 ) 2 )+e-^). (2.8) 

Hence, the new field variables play the role of a normal fields' variables for the problem 
at hand. In these variables the metric acquires the form 

ds 2 =e^dt 2 -e- 2 ^ 2 - 2 ^dr 2 -p 2 e- 2 ^(de 2 +sm 2 ed ( p 2 ) =e 2 *dt 2 -e 2 *dl 2 (2.9) 

where ip(r), <p2{r) arid (p(r) are still unknown functions of the variable r. 

Obviously, the variable ip describes the Fock conformal transformation to the 3D space 
with infinitesimal distance dl F . The variable (p is not a dynamical one and fixes the radial 
gauge. We define a basic radial gauge (BRG) via the relation <*Pbrg(t) = 0. In BRG the 
coefficients of the diagonal kinetic term in (12. 8p are constant. 



2.3.3 Solution of Einstein Equations for Single Massive Point Source 
a) Distributional form of EE for Point Particle. 

Let us consider EE (II. 7p . rewritten in the form 

RZ-S7t(t<?-±T6^ =0. (2.10) 

Since the energy-momentum tensor (I2.2p of the problem is a distribution: T£{x) G 
D''M ( - 1 ' 3 ^{g^ u (x)}, the correct mathematical treatment requires to consider tensor- valued 
test functions ^/^(x) G r DM.^ 1 ' 3 ^{g fMU (x)} [20] and to rewrite equations (12.101) in the form 
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Jm 



(1,3) 



d'x^MxJl - Svr - \TS?j\ %{x) = 0. (2.11) 



For a static problem the expression in the lhs in equation (12.101) does not depend on 
the variable x° and we can use test functions of the form ^(x) = ^(r)x^(a; ), where 
J_2 d x°xfa°) = 1 and ^(r) G DM^^r)}. Then equation (ETTTl) reduces to 

J dh^g tt {v)\ 3 g{v)\ (jtf - 8tt (t? - ^Tsty ^(r) = 0. (2.12) 

These equations must be fulfilled for any test unctions ip^(r). 

b) Solution of the relativistic constraint equation (RCE). 

Taking into account: 

i) the only nonzero component of the energy momentum tensor (12.21) Tq = M5 g (r); 



ii) the expression R$ = -^=/\ g (^fgTt), which is valid in the static case, where A s 

Jj=^i (\/l 3 5 , l5 ,l '' ( ^i) ^ s ^ ne Laplacean in M^jg^r)}; and 

iii) using the specific test functions of the type V^( r ) = i[>(r) 5q5®; 



we obtain the RCE on the initial conditions in the form 



/ d 3 ry/\ 3 g(r)\A g (V&t) V(r) = 47rm^(0). (2.13) 
Here emerges a new constant 

m := MVHO). (2.14) 

In spherical coordinates one easily finds ^/| 3 p(r)| = a/— gw ( r ) p 2 (r)/r 2 . Then using 
the normal field's variables (12.51) and 3D Euclidean space notations, we reach the final 
form of the RCE: 

p 2 J d 3 r V- (Ve^V y~JJ ^( r ) = 4 W>(0). ( 2 - 15 ) 

Its solution determines the dependence of the function (p(r) on the radial gauge function 
(p{r) in the form: 

rr 

(p(r) = ™ dre-^ r \ (2.16) 

P Jr x 

The value > of the radial variable, used in formula (12.161) . defines the place, where 
V ? ( r oo) = 0, i.e., where gtt( r oo) = 1- This value obviously may depend on the choice of 
the gauge function (p(r). The value of the <p(r) at the place of the point source is 

^(0) = -— / dre~^ r) < 0. (2.17) 
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Now we see that the solution of RCE translates the differential 3-form uo 3 : = 
d 3 r y/\ 3 g(r) \A g (y/g~tt) on M^{py(r)} to the distribution- valued differential 3-form uj 2 := 
d 3 v A (— ) = d 3 rm§(v), defined on the Euclidean space E^{5y}, i.e., on the solution of 
RCE we have the relation 

u 3 = u 3 . (2.18) 

One can consider this relation (12.181) as a new form of the RCE. 

The correspondence between the spaces M^{^y(r)} and ~E^{Sij} was stressed in 
Section 2.1. The extension of this correspondence, obtained here, is the geometrical basis 
for application of the mathematical theory of distributions in the massive point particle 
problems in GR. 

c) Solution of the Rest of EE. 

Since the other components of the energy-momentum tensor of point particle are zero 
in its proper frame, one can use for them any of the forms (11.71) . I2.10p . or (12.111) of EE, 
and we arrive at the following ordinary differential equations for normal field variables: 

^ + ^3 = 4^"^ (2.19a) 
P 

(^f - (^f + \e 2 ^ 2 -^ = (2.19b) 
P 



Using the relation (12.161) . we can exclude the gauge function p> from this system, 
its radial-gauge-invariant form: 

dp> 2 m 2 
\dp J m z 

Note that meanwhile we have excluded the radial-gauge-dependent variable r replacing it 
with the radial-gauge-independent one - <p, which now plays the role of the independent 
"radial" variable. 

The first equation (I2.20al) can be integrate immediately in quadratures. The second 
equation (I2.20b|) imposes a constraint on the two integration constants in the general 
solution of the first one. Thus we remain with only one integration constant p,^ in the 
solution of the system (I2.20[) : 

(p 2 ((p) = In ( . ry— r 

VI smh(v9 - v?oo) 



mexp(— (p) 
sinh Up - poo 



(2.21) 



d) Fixing the emerging constants in the general solution of the problem. 

The second expression in equations fl 2 . 2 1 [) is derived using formulas (12.51) . It shows that 
p(v 9 oo) = oo. Hence, the value p^ corresponds to the physical infinity, where the space- 
time is asymptotically flat and we must have g tt = 1- Thus we see that the value p^ 
must be reached for the value of the radial variable r^, i.e. the relation p^ = pir^) = 
must take place. For the value of the radial variable r ro the space- time M^ig^ix)} 
is asymptotically flat. This value corresponds to the physical infinity. Hence, under our 
conventions, described in the Introduction, the physically admissible interval of the values 
of the radial variable is r e [0,^]. 

As a result we remain with only two arbitrary constants m and in the solution of 
the whole system of EE for the massive point particle. The 4-D metric in normal field's 
variables acquires the final gauge invariant form in which its coefficients are functions 
only of p> in the role of a radial variable: 

,9 9,„, 9 9,n ?f dp 2 d6 2 + sm(6) 2 dd) 2 \ 

ds * = e ^ dt 2 - m 2 — -f — + . \ 2 V ■ (2.22) 

\smh(<yj) 4 smh^py J 

It's remarkable that in the metric (12.221) appears only the integration constant m. As a 
result only the value of this constant will influence the local dynamics of any test particles 
and fields, which probe the metric in the space-time of single massive point source. This 
important conclusion is independent of the choice of radial variable r, i.e. it is gauge 
invariant, as well as the whole equation f!2.22[) . 

According to the definition (12. 14j) of the constant m and the formulas (12. 5p we have 
p>{0) = \n(m/M). The relation (I2.17P shows that <^(0) < 0. Thus we obtain, that in the 
GR massive point particle problem the variable p varies in the interval p> G [ln(m/M), 0]. 
After all, the bare mechanical mass M shows up in GR massive point particle solution, 
defying the interval of the physical values of the variable <p. 

It is convenient to introduce the mass ratio 



q = m/M e (0,1). 



(2.23) 



Then in the problem at hand the basic quantity ^fg~tt = exp(</?) varies in the physical 
interval ^/gj t G [g, 1]. 
The final formula 

. , 2m 2m 

P(<f) = 7^ > 1 2 > 2m 2 - 24 

1 — exp(2{p) 1 — g 2 

shows that the luminosity variable p in the gravitational field of massive point particle 
cannot take values, less than > pa, since 2m = pc, as we shall see in the next 
subsection. This is in strong contrast to the situation with Schwarzschild solution in 
Hilbert gauge (jl.6p and in complete accord with Dirac's suggestion [2~T] . 

d) Fixing of the physical radial variable. 

The previous consideration gives a correct mathematical ground for our conclusion 
about the topology of the space-time of single massive point particle problem in GR, 
as described in Section 2.1. Indeed, it is easy to obtain from EE (I2.10p the 4D scalar 
curvature of space-time with single massive point source: 



and its 3D scalar curvature 



4 R= -8vrM4 3) (r), (2.25) 



3 R = 167rM^ 3) (r). (2.26) 



Clearly, the last equation is the concrete form of the RCE (11.151) in the point particle 
case. Hence, the 3D space M^^r)} has a strong singularity at the geometrical point 
r = 0, where the massive matter point is placed, and the solution of EE cannot be extended 
behind this point, both from physical and from mathematical reasons. The space-time 
of the problem M^ 1 ' 3 ^{g fJjl ,(x)} has a singular line - the world line W of the massive 
matter point. This singular line must be removed from the manifold M^ 1,3 ) {g^ ll/ (x)} and 
we remain with the topology, described in Section 2.1. 

The simple correspondence between the Riemannian space-time of point particle and 
Minkowskian space-time gives a good reason to adopt, as much as possible, the basic 
properties of the Minkowskian radial variable r for the Riemannian CclSG, clS described 
in Section 1.2.2. Then, according to results in the previous Section, we have to fix the 
gauge function tp(r) in such a manner, that as a result of relation (12.161) we will have the 
mapping: 

[ln(m/M),0]^^![0,oo] r (2.27) 

A very important additional requisite of the mapping (I2.27P is the requirement to 
preserve the number and the character of the original singularities of the solution (I2.22p 
in the whole compactified complex domain of variable <p. The mapping (12.271) is allowed 
only to change the positions of these singularities in the compactified complex domain of 
variable r. 

The only way to fulfill this requirement is to use a fractional- linear function ip{r) = 
with some constant coefficients a, b, c, which are unambiguously fixed by the mapping 
(I2.27P in the form: ip(r) = r+b / lng - Taking into account that at r — ► oo the asymptotic of 
the function (p(r) ~ b/r yields an asymptotic g t t ~ 1 + 2b/r, one sees that the standard 



comparison with the real observations imply b = —m Kepler, where mxepier is the Keplerian 
mass of the particle, as observed by a distant observer. Then the formula (I2.24[) gives 
lim = m/rriKepier and the Edington's coherent scale condition: lim = 1 (see 

Section 1.2.2) fixes the value of our integration constant m = rriKepier- Hence, in this 
physical gauge the function 

<p(r) = - , ,™ , = - , M f^ p = <p{r; M, R, g) (2.28) 

|r — r \+m/m± |r — r | + R 

presents a proper GR generalization of the Newton potential (11.21) if Newton (r) = — 

of matter point with Keplerian mass m, placed at the position r , which describes the 

fundamental solution in Newton theory of gravity^. 

Now we see that the ratio g = m/M G (0,1) describes the relativistic gravitational 
defect of mass for massive point particle. It was introduced for the first time in [16] , where 
the solution of the problem was derived using a different mathematical technic. As a final 
result we obtain a two parameter family of solutions (I2.28|) to the massive point particle 
problem in GR. This family can be parameterized by any two of the three constants m, 
M, and g. 

In the formula (12.281) we use a short notation R = mj In ~ = M j y~ In ~J for the GR 
correction to the Newton potential. Further on we shall refer to the correction R as "a 
relativistic shift" in the Newtonian potential. The potential ip Nawton (r) can be derived as 
a limit R — > of the relativistic one Therefore we shall refer to the very function 

<p(r) (I2.28P as "a relativistic gravitational potential". From mathematical point of view 
the relativistic gravitational potential defines the Fock conformal mapping. It is clear 
that this potential plays a basic physical role in the relativistic theory of gravity. 

2.3.4 Some Remarks on the Non-Relativistic Limit c — > oo. 

One can expect that in the non-relativistic limit c — - > oo our GR solution for single particle 
will reproduce the results of the Newton theory. For study this limit it is necessary to 
restore the physical units in corresponding formulae. Then we obtain for the mass ratio: 

g = e^ 0)/c \ (2.29) 

and 

QNewton M e ^(0)/c 2 

^ = ~|r - r | - G«»»Ie^«)/ c 7^(0) (2 ' 30) 
- for the relativistic potential. 



lr The form of the Newton potential can be derived, following the same consideration with only one 
difference: the physical values of the Newton potential cover the whole semi-constrained interval (y9 Nc "' t °™ g 
(—00,0]. This interval has to be mapped onto the interval [0, oo) 3 r by fractional-linear function 
<p ewi ""* (r) and this gives (/9 Nc ™ t °" ( r ) = 

2 One easily obtains the following instructive estimates: a) If g < l/\/e ~ 0.60653, then R < pa- b) If, 
according to Birkhoff 's theorem, one applies the formula (|2.28[) outside the spherically symmetric body of 
finite radius rg, the quantity R will not exceed the radius of the body rg, when qb < ex P(~PGl r B)- This 
restriction is very weak, since Pg/i"b < SC 1 for real bodies, and for any of them 1 — exp(— pg/ t b) <^ 1- 



Unfortunately, at present we do not have a theory of the relativistic collapse, which 
has to describe in detail the origin of the relativistic gravitational defect of mass of a 
single point particle and the value of the mass ratio g. Here we are considering g just as 
an additional free parameter of the class of point particle solutions, studied in the present 
article. If one considers, instead, as a free parameter ip(0) and assumes that it does not 
depend on the velocity of light, in the limit c — > oo one obviously obtains from the relation 
(TQj]) the limit: 

lim^^ g = 1. (2.31) 

This result sounds physically right. Indeed, one expects that in the non-relativistic limit 
c — > oo the gravitational mass defect will disappear and we will return back to the Newton 
theory of gravity with m = M. 

Despite of this physically reasonable result, the assumption that y?(0) does 
not depend on the velocity of light c gives a wrong limit /im c ^oo (p(r) = 
_QNewton M j _ ro | _ QNewton M/(p(0)) in the formula f[2T30|) . 

To obtain the physically right results in the both cases, one has to assume that: 

i) actually the quantity <^(0) depends on the velocity of light c in some specific way, 
and 

ii) the unknown at present function ip(0, c) fulfills simultaneously two additional con- 
ditions: 

lim^ <p(0, c) = <p Newton (0) = -oo, 

Zzm™ (<^(0,c)/c 2 ) = 0. (2.32) 

As a result of these conditions, which are obviously compatible, one obtains both the 
relation ( 12.311) and the right non-relativistic limit 

s~iN 'ewton tut 

Um^ <p(r) = --. = ip Newton (r). (2.33) 

|r - r | 

One can hope that the future theory of the relativistic gravitational collapse, accom- 
panied by a proper treatment of the gravitational mass defect, or some other additional 
considerations, will be able to derive the precise form of the function (p(0, c) and to con- 
firm the physically natural relations (12.321) . In the present article we will assume these 
relations to be fulfilled. 



2.3.5 Some Remarks on the Properties of Static Fundamental Solutions 

1. On the three dimensional form of the fundamental solutions of EE. 

An unexpected and remarkable feature of the relativistic gravitational potential (I2.28j) 
is that it has a finite negative value <p(r ) = In g at the place of the very point source r . 
This unique property is in a sharp contrast to the case of the Newton potential ^ Newton ( r ) 
(11.21) . which diverges as — m/|r — r |, when r — > r . Thus we see that in GR we have a 
self-regularizing mechanism for gravitational interaction, based on the influence of matter 
on the space-time curvature. In the article [22J we have shown that the same phenomenon 
comes into being in GR electrostatic problem of single massive point charge. 



As a result of the relativistic self-regularization all components of the metric tensor 
y/gtt — \/—g r r = s ^h(ip) 2 ' P = sinh(-y) * n spherical coordinates are regular at the 
place of the point source r : 

/ — i — 7 / 7 — r 4/f)ln^\ 2 2m , nn A\ 

V9tt{ro) = Q, V-9rr{r ) = - I J , PM = ~r~Z 2 ' ( 2 - 34 ) 

Now it becomes clear that the singular term #( 3 )(r — ro) in the lhs of Eq. (I2.15P 
originates from the singularity of the 3D Cartesian determinant: 

a/IM 1- )! = V-9rr(r)p(r) 2 /\r - r | 2 ~ l/|r - r | 2 , when r ^ r . (2.35) 

The singularity of metric coefficients at the place of the point source of gravity becomes 
transparent in Cartesian coordinates. Indeed, one can write down the 3D distance in the 
tensorial form dl 2 = —dr 3 g(r) dr, using the 3D Cartesian metric tensor 

— 3 g(r) = — (I — e r ® e r ) — g rr (r) e r <g> e r , where e r := -. (2.36) 

|r - r | 2 |r- r | 

The components of this tensor are obviously singular at the point r = ro- This is precisely 
because in the specific geometry, defined by fundamental solutions of EE, we have p(r ) > 
0. The expression (I2.35P defines the square root of the determinant of tensor (12.361) . 

2. One-dimensional-like representation of the fundamental solutions of EE. 

The Dirac 5-function is a linear functional. Its representation depends on the class of 
the test functions in use. We can take advantage of spherical coordinates in description of 
the test functions of the GR point particle problem. Starting with Cartesian coordinate 
test functions G D / {M ( - 3 - ) }, in spherical coordinates we obtain a specific class of test 
functions ip sp h(r, 9, 0) := ip(re r (9,(j))) G D' h {Mi^ + x SO (3)}. These must be distinguish 
form the arbitrary test functions ifj(r,6,4>) = i/)(r,e r (9,(/))) G "D'{M.^ + x SO(3)} on the 
manifold M( 3 ) = x SO (3). Now r G is considered as an independent variable, 

not just as a short notation for |r| := a/x 2 + y 2 + z 2 . The difference between the functions 
Tpsphir, 9, 0) : = VK re r(#, 4>)) an d i[>(r, 9, 0) = if)(r, e r (9, 0)) is of critical importance for our 
problem, since ipsph{ r — 0,0,0) := ip{0) = const for any values of the angle variables 9 
and 0. The functions -0( r , 9,4>) = ijj(r,e r (9,<j))) do not have such property. Instead, the 
functions ip{r = 0,0,0) = -0(0, e r (9, 0)) ^ const keep the dependence on the variables 9 
and 0. 

Let us use the class of test function D' h {M.^ + x 5*0(3)}. The standard restriction of 
the Euclidean Laplacean A(p(r) = - d 2 2 (np(r)^ on the functions, which depend only on 
variable r, brings us to ID formulation of the problem. One can write down the solution 
of the Eq. (12.151) with added point source of gravitational field, described by function 
S(r), in the following one- dimensional form (see for details [16J): 

^) = i^(i-^e(^)) 

Here we are using the Heaviside steep function B(r) with regularization O(0) = 1. 

The fundamental solutions of EE were found for first time in [16] in this form. Its 
advantage is that it makes transparent the jump in the derivatives of the metric coefficients 



(2.37) 



in spherical coordinates. This jump reproduces via the Einstein tensor the ^-function in 
the rhs of EE with point source. 

In the case of ID representation (I2.37P the form of the metric (12.221) must be considered 
as valid only in the vacuum domain, outside the point source. In this domain the form 
(I2.22p does not make difference between 3D and ID representation of the fundamental 
solutions. 

2.3.6 On the Choice of Radial Gauge in the Single Particle Problem 

The above consideration solves on a clear theoretical basis the longstanding problem of 
the choice of radial variable r for point source of gravity in GR. The unambiguously 
obtained physical radial variable r is obviously a preferable one, both from mathematical 
and from physical point of view. 

The singularities of the metric coefficients in the whole compactified complex plain 
are placed at the positions cp n = inn, n G Z. The points of finite n are poles and the 
infinite point \n\ = oo is an essentially singular one. 

The singular points of the solution in the whole compactified complex plain C r of the 
physical variable r are of two essentially different types: 

1. The place of the point source of gravity at r = where the curvature of space- 
time has a strong singularity, proportional to 5^ 3 ^(r). This singularity is seen in the 
differential 3-forms (12.181) . too. Surprisingly, the relativistic potential f!2.28|) and the 
metric coefficients, when written in 3D form, are regular at this point. 

2. In Einstein theory of gravity an unavoidable singular points of the metric coefficient 
gtt{ r ) are the (complex) points r n = —R + i^, n G Z. These singular points are placed in 
the nonphysical domain of the physical variable r. For finite n G Z the singular points are 
poles of the 3D metric coefficients in the expression (12.221) . The sequence of these singular 
points has a limiting point r = —R for \n\ = oo. This is a real essentially singular point, 
which is not an isolated one. In contrast, the real singular point r = —R is a simple 
pole of the relativistic potential (12.281) . which has no other singular points. This pole is 
placed in non-physical domain r < and corresponds to the real pole at r = of the 
gravitational potential in Newton theory of gravity. 

The multiplication of this single simple pole of the relativistic potential ip(r; M, R, g) 
to an infinite series of singularities of GR metric coefficients is produced by the specific 
exponential mapping (12.221) of <p(r; M, R, g) onto these coefficients. This is a specific 
feature of the relativistic description of gravitational field of massive point particle and 
may be considered as a price, one has to pay, for the self-regularization mechanism, 
discussed in the previous subsection. 

Under gage transformation to some other radial variable r ot h eri related with the phys- 
ical one r by a coordinate transformation r other = r other {r) of general type, which is not 
fractional-linear one, in the solution of the problem will appear additional nonphysical 
singularities in the corresponding compactified complex plain C rother . 

From the relation (12.161) one easily obtains the physical gauge function ^ p hys(r) = 
21n((r + R)/p), which is compatible with boundary conditions of the problem and with 
all additional requirements on the physical radial variable r, as formulated in the Intro- 
duction. 

We shall call a regular gauges of the problem all gauges, for which the integral (I2.16P 
makes sense in the physical interval r G [0, Too]- This is just the necessary and sufficient 



condition for the radial gauge to be compatible with the boundary conditions, coded in 
the Dirac 5-function in the RCE of the massive point particle problem. 

This condition fixes a large class of admissible gauges for this problem. One of them 
is BRG (p(r) = in which the relativistic potential has the form brg(t^) = m ( r — 
^oo)/p 2 , r e [0,7-00]. 

Between the regular gauges for the one-particle problem in GR is the gauge by Droste 
[I]. It has a clear geometrical meaning, since in this gauge the radial variable r measures 
the radial geometrical distance in the 3D Schwarzschild metric. This gauge reproduces 
only a very special value of the mass defect ratio g = (^5 - l)/2 « .6180 [16J. Quite 
curiously, under such geometrical choice of the radial gauge g equals precisely the famous 
mathematical golden ratio. 

All other known radial gauges, probed for spherically symmetrical static solutions of 
EE and described in [16], are not regular. Therefore they cannot be used for solution 
of the point mass problem. As we have seen in the previous subsections, only a certain 
combination of gauge function (p(r) and corresponding form of the relativistic potential 
f(r) can obey the specific boundary conditions of this problem. 

For example, in the most popular at present Hilbert gauge: p H (r) = r the static spher- 
ically symmetric problem has a relativistic potential ip H (p) = In a/1 — Pg/p and a radial 
gauge function (p H (p) = In — Pg)/(p) 2 )- Hence, for Hilbert gauge the integral (I2.16P 
diverges logarithmically: 1) when p approaches the center p = 0, where the point source 
of gravity must be placed in this gauge, if one insist on the point particle interpretation 
of the this form of the Schwarzschild solution; 2) when p approaches the event horizon 
p = pq. In addition, the value of this integral becomes an imaginary number for p < pc- 

This means that Hilbert gauge is incompatible with the specific boundary conditions 
for EE in presence of massive point particle. Therefore one cannot use the Hilbert gauge 
to solve the point particle problem in GR. This gauge yields the well known nonphys- 
ical singularity at the point p — 0, i.e. on the boundary of the to-be-physical domain 
of the radial variable p e [0, 00). More over, the meaning of the variable p radically 
changes in the interval [0, po\- Here it plays the role of a specific time variable and the 
point p = describes the future infinity of the internal time U n = x — 1/x G (—00, 00), 
where x = p+pc In (Ip/Pg — 1 1) is the Regge- Wheeler "tortoise" coordinate in the interior 
of Hilbert solution [23]. It becomes clear that even if we will be able to reproduce math- 
ematically a term ~ 5{p) in the rhs of EE (see the articles by P. Parker, by H. Belasin 
& H. Nachbagauer, and by J. M. Heinzle & R. Steinbauer in [15J), its interpretation as a 
source of the gravitational field and the curvature of the Schwarzschild solution is physi- 
cally unacceptable. Such term may describe only a 5-shaped - with respect to the time, 
"impulse" at the time instant p = (=>■ t in = 00) and has a complete unclear physical 
meaning. In any case it is not able to describe the usual physical 3D-space-point source 
of static gravitational field. 

2.4 The Total Energy of the Aggregate of Massive Point Source 
and its Gravitational Field 

In the problem at hand we have an extreme example of an "island universe". In it a 
privileged reference system and a well defined global time exist. It is well known that 
under these conditions the energy of the gravitational field can be defined unambiguously 
PQ. Moreover, we can calculate the total energy of the aggregate of a mechanical particle 



and its gravitational field in a canonical way, considering the corresponding ID varia- 
tional problem for total action (12.11) in the spherically symmetric static case [16]. The 
canonical procedure produces a total Hamilton density "K t ot = S a=1)2;A1 =t /r Va,^ — ^tot = 
| [—p 2 ip' 2 + p 2 ^2 — e 2v2 ) + Me v 5(r). Using the equations (12.191) . one immediately ob- 
tains for the total energy of the GR universe with one point particle in it: 



This result completely agrees with the strong equivalence principle of GR. The energy 
of the static longitudinal gravitational field, created by a point particle at rest is a negative 
quantity: Eqr = E tot — E M = m — M = — M(l — g) < 0. Since both matter point and its 
gravitational field have nonzero proper energies, this result proves that the ratio g must 
belong to the open interval (0, 1), see [16] for more details. 

The above consideration gives a clear physical explanation of the gravitational mass 
defect of a point particle. 

3 Quasi-Linear Superposition Principle for Static 
Fundamental Solutions of Einstein Equations 

3.1 Justification of the Quasi-Linear Superposition Principle for 
the Relativistic Gravitational Potential ip 

3.1.1 Some General Arguments 

As a result of the GR dynamics, the relativistic constraint equation (RCE) is a restriction, 
which will be fulfilled at any time instant t, if it is valid at the initial time instant t = 0. 
Therefore it is enough to solve RCE only at the initial time instant t = 0. 

One can use this specific feature of the RCE to simplify it, imposing special additional 
conditions at the initial time instant, which can not be fulfilled during the further evolution 
of the physical system. Thus one can impose different physical conditions on the initial 
state of the system under consideration. 

For example, as stressed already by Misner in [17] and a bit later by Fock in pQ, it 
is impossible to have a permanent static solution for N-particle system in GR, if N > 1. 
Nevertheless, it is possible to find an initial-instant-static solutions of the problem with 
any number of particles N. As a rule, the initial conditions will contain some initial 
amount of gravitational waves. One of the basic open problems is how to exclude the 
presence of initial gravitational waves. One can expect that for such solutions RCE will 
take its simplest form. 

It turns out that under the conditions d t g^ v \ t= Q = on the time-derivatives of the 
metric and some weak additional constraint on d 2 2 gfiu\t=o the RCE, together with other 
EE, yields a quasi-linear equation for the relativistic potential of any mass distribution. 
This equation can be considered as a relativistic analog to the linear equation (jl.ip for 
Newton potential in the classical theory of gravity [T7] . 

Let us consider a system of N point particles with bare masses M^, A = 1,...N 
at positions r^t). Here we suppose to work with space-time manifold, which allows an 




■oo 



"Ktotdr = m = gM < M . 



(2.38) 



existence of global time t. The energy-momentum tensor of such system at time instant 
t is: 



N 

TZ(t,r) = J2 M A 3) (*-*A(t))<(t)u Alx (t). (3.1) 

A=l 

According to the relativistic nonlinear superposition principle, described in Section 
1.2.4, the metric g^y{t,r) of this N-particle problem is given by formula (11.131) . This 
formula does not yield any simple practical results. It describes in a formal way the 
solution of the very complicated GR problem under consideration and demonstrates its 
existence and uniqueness under proper boundary conditions. 

As we shall show in this Section, in contrast, one can introduce a simple quasi- linear 
superposition principle for the relativistic gravitational potential <p(r). 

We shall consider the special case in which at the initial time instant t = all particles 
are at rest and have initial positions r A (0) = r A . Then 

T;(0,r) = (f>^f)(r - r A ) \ S^l. (3.2) 

According to articles [T7j, the generalization of the covariant form of RCE (12.261) for 
N point particles at rest and under additional conditions dtg^ u \t=o — 0, i.e. in the case of 
the N-particle instant- gravistatics, is: 

JV 

3 R = 16vr^M A ^ 3) (r-r A ). (3.3) 

A=l 

The proper generalization of the relation (12.131) . which follows in the instant- gravistatic 
case from the EE for N-particles reads: 

r N 

/ d 3 r y/\3g(r) \A g (^) ^(r) = 4vr V m A ^(r A ). (3.4) 

Jm(3) 

Here m A = M A ^J g tt (Y A ) are the corresponding Keplerian masses of the point particles. 

A remarkable feature of the equation (13.41) is that under conformal Fock transforma- 
tion: 

g l3 = e-^h l3) ^' = e 2 ^, ^g~\ = e -^yfh\, (3.5) 

one obtains the quasi-linear equation for the relativistic potential of the N-particle prob- 
lem: 

N 

di (yfh\h ij dj^ =4ttJ2 m a e VA 5 {3 \r - r A ), where <p A := <p(r A ). (3.6) 

A=l 

Here are functions which define the unknown metric in the Fock conformal space of 
the N-particle case. In Eq. (13. 6p we have used the substitution ^/gTt = e^. 

After the pioneering work by Lichnerowiz, one usually supposes the metric hij to be 
conformally flat (see in |T7]). This leads to a well known superposition principle in the 
case of instant gravistatics of N Schwarzschild black holes |17j . 



Today it is well known that the conjecture of conformal flatness of hy is too restrictive 
and does not allow one to obtain the solutions of real physical problems (see the review 
article by Cook in [17] and the references therein). Unfortunately, at present we do not 
know an alternative assumption, which fixes the metric hy in a physically acceptable way 
for the case of N matter bodies. 



3.1.2 The Superposition Principle for the Potential ip 

Here we outline a complete different approach to the problem at hand. It is based on 
a specific superposition principle in GR. This new approach may turn to be a more 
physical alternative to the conformally flat one, mention in the previous Section. The 
novel superposition principle may play the role of the additional physical requirement, 
needed to select the proper initial conditions for the GR N-body-problem between all 
mathematically admissible and formal initial conditions. 

Having in mind the inhomogeneous quasi-linear equation ( 13. 6ft . it seems natural to 
define its solution, we are looking for, by the formula 



N 



y(r;n,...,r N ) = -^ (3-7) 



A=l 



|r -r A \+R 



It follows a simple quasi-linear superposition principle for the static relativistic potential ip. 
Here the Keplerian masses are uia = QaMa and the mass-defect ratios qa(*i, ■ ■ ■ , fn) = 
exp <pa{ta] ri 3 • • • , i"n) of the A-th particle in presence of the other massive points can be 
obtained from the self-consistency condition for the relativistic potential 



N 



ip(rA',r 1 ,...,r N ) = ^2c B (p(\r A -r B \,M B ,R B ), (3.8) 

B=l 

V (\r A -r B \,M B ,R B ) = ^— - . (3.9) 

\*A — *B\ + J^B 

Note that the procedure, based on the relations (13.71) and (13.81) . represents the only 
way to construct a relativistic gravitational potential ip(r; r±, . . . , r^), (and corresponding 
metric) of a system of N point particles with the following 

Fundamental property: For a system of N particles at finite and nonzero mutual 
distances the total variety of singularities of the relativistic potential in the whole com- 
plex domain of the variables is just a superposition of the singularities of the relativistic 
potentials of the separate matter constituents of the system. 

In other words, joining several point particles in a gravitationally interacting system, 
we remain just with the singularities of all independent particles. As a result of the 
integration of the independent particles in a joint interacting system no additional new 
singularities emerge, as well as no old singularities disappear in the whole complex domain 
of space variables. In addition, the singularities of the separate particles do not change 
their character. 

From analytical point of view one can reach such simple preservation of the singularities 
only by constructing a linear combination of the corresponding analytical functions with 
some constant coefficients. The uniqueness of this construction, up to the choice of the 
constant coefficients, is guaranteed by the corresponding theorem of the complex analysis, 
which stays that every analytical function is unambiguously defined by its singularities. 



In the quasi-linear superposition principle (13. 7\\ we are using the single point particle 
solutions (12.281) . denoting by = exp {<p°£) the value of the mass defect ratio of the 
A-th particle in the case tab '■= \ r A — y b\ — * °o for all B ^ A, i.e., in the previously 
considered case of a single massive point particle in the whole universe. The constants 
Ca = Qa/ ' Q°a > are unknown and have to be justified. 

The above consistency condition ( 13. 8ft yields the following basic nonlinear algebraic 
system of N equations: 



for the mass defect ratios qa = exp (<£>a) of the A-th particle as a member of the N-particle 
system in finite space range, i.e. when all tab < °o- The system (13.101) is regular one 
when Ra ^ for all values of A = 0, N. 

In the N relations (I3.10p we have too many free parameters, which have to be fixed 
using some proper physical assumptions. Taking into account that the bare mass Ma and 
the relativistic shift Ra = Ma/ (jpz In -^j are inner characteristics of the very A-th point 

particle, we can suppose both of them to be constants, whose values are independent of 
the N-particle configuration. 

Then we see that our quasi-linear superposition principle (13.71) leads, after all, to a 
definite mathematical formulation of a novel relativistic Mach-like principle. It states that 
the Keplerian masses rriA = wi^r^, . . . , rjv_i,jv) of the bodies depend on the mass distri- 
bution in the universe, in contrast to their bare masses Ma, which remain independent 
of matter distribution. These essentially different properties of the masses rriA and Ma 
seem to be natural in the relativistic theory of gravity and may have important physical 
consequences. 

The independence of the bare particle masses Ma of the system configuration seems to 
be quite natural assumption. In addition we will assume the bare mass to be an additive 
quantity, i.e., the total bare mass M of the composite system of particles is just the sum 
of the bare masses of the constituent particles: M = J2a=i Ma- 

More speculative, from physical point of view, is the requirement of the independence 
of the relativistic shifts Ra from the particle configuration. In the present article we test 
this assumption as a way to restrict the number of the free parameters in the problem at 
hand and study some of its consequences, leaving for future developments its justification 
on a more profound physical basis. 

3.1.3 The Non-Relativistic Limit of the N-Particle Potential ip 

The restoration of the correct physical units draws an additional light on the proposed 
quasi-linear superposition principle (13. 7p . In physical dimension-full quantities Eq. (13. 7p 
acquires the form 




(3.10) 




(3.11) 



and the consistency conditions (13. 10p read: 



N /^Newton Jl r in B 

r , = e(—?)/=V? " £ " ■ ■ A=1,...,JV. (3.12) 

BjtA r AB -G Newton M B e^/ c /cp% 



The last formulas show directly that under the two assumptions (12.321) in the limit c — > oo 
one obtains precisely the Newtonian superposition principle (11.31) . as it should be for the 
correct relativistic generalization of the non-relativistic theory. This observation increases 
our confidence in the approach, based on the superposition principle (13. 7L 
To proof this statement one should take into account that: 

a) lim^oo <p°g = — oo and lim^oo (<^/c 2 ) = - according to the relations (12.321) . 

b) Using the last relations in the consistency condition ( 13.121) one obtains easily first 

lim^oc ((Pa/c 2 ) = 0, (3.13) 

and then 

lirric^ootpA = -oo. (3.14) 

c) As a result qa = e lpA ^ c2 — > 1 and niA = QaMa — > Ma when c — > oo. 

3.1.4 The Case of Continuous Mass Distribution 

The generalization of the equations (13.71) for continuous distribution of mass is straight- 
forward and reads: 

<p(r) = - / rfV , ^ , ■ - - ■ (3.15) 



M( 3 ) 



|r - r'\ + R(r') 



Here /i(r) is the density of bare mass M. In the case of continuous mass distribution the 
relativistic shift R is R(r) = /i(r)/x°°(r). The local density x°°( r ) °f the quantity x°° is 
the second independent function, needed for description of the density of the relativistic 
point potential in the case of continuous mass distributions. 

The relation ( 13.151) generalizes and replaces the non-relativistic superposition principle 
(11.41) for continuously distributed masses. The last can be derived from equation ( 13.151) . 
taking the limit c — ► oo precisely in the same way, as in the discrete case. 

From mathematical point of view the relation (13.151) is a nonlinear and nonsingular 
integral equation for the relativistic potential <p(r): 

ip{r)= I d 3 v'K(r,v')e^ r,) (3.16) 

with a nonsingular kernel 

For continuous distribution of identical particles with fixed x°° — const one has to put 
R(r) = /i(r)/x°° in (I3.17p . thus remaining with only one given function /x(r) > in the 
kernel -^(r, r'). The function fi(r) reduces to a given constant fi = const for homogeneous 
mass distributions. Hence, in the last case R = const, too. 



3.2 Some Basic Solutions of the Mass Defect Equation 



In this Subsection we are testing our basic assumptions, described in the previous Sections, 
considering both the cases of a few point particles, and of continuous mass distributions. 
For this purpose is convenient to introduce the quantity 



M 11 
x :— — = - in - > 0, 
R g g 



(3.18) 



which turns to play a basic role in our considerations and demonstrates some simple 
properties. Making use of the Lambert function W(z) [21] , i.e., the solution of the equation 

We w = z W = W(z), 



one obtains, solving the equation (13.181) with respect to the mass ratio g, the basic relation 
(see Fig. HD 

e = g( x ) = e~ w{x \ (3.19) 
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Figure 1: The dependence (13.191) of the mass ratio g on the variable x. 



Now one can write down the basic equations (13. lOf) in the form: 

N 



In 



Xb Qb 



\ 1 + r AB /R B ' 



(3.20) 



We shall call this equations mass defect equations (MDE). One immediately obtains from 
MDE two important consequences: 

3.2.1 A Few Particle Solutions 

1. Let us consider first the two-particle solutions (N=2) of the MDE: 
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In 






\Qi 






In 






V^2 



1 + r 12 /R 2 



eft 



1 + r 12 /Ri 



+ x^g 2 



(3.21) 



For them we have two interesting limiting cases: 

i) Total decay of the two-particle system, when r 12 — > oo. Then obviously 

Qi(r 12 ) -> §T, ^2(^12) -»• Q?, 

as one expects. 

ii) Merger of two particles r 12 — > 0: 

eiu2(ef > , 2 °°) = exp ( - W (xf + x?) ) . 
Note that for the quantities x we obtain a simple linear superposition: 

x lu2 = x™ + x?. (3.22) 

As a result 

1 = £ + £=:/4\, ^^^[0.1], (3.23) 



#iu2 #2 i?i W/ ' M x + M 2 



Introducing the quantities 



SM:=frf; and M ;= g^g, (3.24) 
Mi + M 2 R1 + R2 



we obtain finally 

l-5i? 2 



1U2 



Here (-R) ar := ^ Xm=i denotes the arithmetic average. In the present case N = 2. 

Now we can derive easily the following basic properties of the function i?iu2 
R lu2 (Ri,R2;M u M 2 y. 

1) R 1U2 (R ± , R 2 ; Mx, M 2 ) = R lu2 (R 2 , R 1 ;M 1 , M 2 ). 

2) R 1U2 {R u R 2 ; M u M 2 ) = R lu2 (R ± , R 2 ; M 2 , M 1 ). 

3) R lu2 (kR 1 ,kR 2 ;M 1 ,M 2 ) = kR lu4 (R u R 2 ; M u M 2 ), VA; > 0. 
5) J R 1U2 (Ri, R2; fcMi, kM 2 ) = R lu6 (R u R 2 ; M u M 2 ), V k > 0. 

7) i? lu2 (R, R;M 1 ,M 2 )=R,VM 1 , M 2 . 

8) R 1U2 (R u R 2 -M,M) = jgg£, or 

-Riu2 (-R) ar = (-R) 9 2 m • 



Here (R) gm ■— \fR\R 2 denotes the geometric average. 

9) i? lu2 (R 1 ,0;M 1 ,M 2 ) = 0. 

10) i?iu2 (Ri, 00; Mx, M 2 ) = ^ = ^JM, where x lu2 = an. 

iii) In the case of two identical particles: Mi = M 2 = M, R± — R 2 — R, at finite 
distance r\ 2 one easily obtains for £>i(ri 2 ) = £2 (^12) = q{tv2) 

,(r 12 ) = exp( - W (l±^x~) ) (3.26) 
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Figure 2: The dependence 03. 26!) of the mass ratio £>(r 12 ) = m(r 12 )/M of system of the two 
identical particles with different individual mass ratios g°° = 0.1, . . . , 0.9; on the distance 
r 12 between them. The distance r 12 is shown in units of R. 



and 



Qiu2{Q 



exp 




The small variation of the mass ratio of each of the two identical particles in the 
system with the change of the distance between the particles is shown in Fig. [2j As 
seen, at distances ri 2 ^> R the measurable Keplerian mass of each particle is practically 
constant, since the value of q is almost constant. 
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Figure 3: The dependence (I3.27f) of the mass ratio difference A^(a; c 
value of x°° for the two identical particles. 



on the individual 



The formula 

Apix 00 ) = g(0) - g(oo) = e~ w ^ - e ~ w ^ (3.27) 
shows that the variation of of the mass ratio of each particle £>(r 12 ) is not bigger then 



~ 0.1408, when the distance between them varies from zero to infinity: r\i G (0, oo), see 

Fig. m 
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Figure 4: The dependence ( I3.28P of the energy release (in units 2Mc 2 ) on the individual 
mass ratio g°° of the two identical particles. 



The energetic efficiency of the process of gravitational merger of two identical particles 
is described by the quantity: 

^ = 1 - ftuaCe 00 )/^. (3-28) 

Here AEq is the energy release in the gravitational collapse of the pair of the two 
identical point particles from infinite distance to their merger. 

2. For three identical particles at the vortices of equilateral triangle tab = t for all 
A, B = 1, 2, 3; A ^ B one obtains: 

g(r) = exp| -w(j±^^) )• (3.29) 

3. For four identical particles at the vortices of equilateral tetrahedron tab = t for all 
A, B = 1, 2, 4; A ^ B the result is: 

g(r) = exp| -Wt^^x-) ). (3.30) 

3.2.2 Some Basic Properties of the N-Particles Solutions 

Consider now a system of N particles: 

i) If the system of N point particles decays, i.e., when tab — » oo for all B ^ A and we 
remain with only one such particle in the whole universe, then it should be. 

ii) In the opposite case, when all tab 0, i.e. when the system N particle collapses and 
they merge into one composite particle, from MDE one obtains Qa = Qb = const = Qiu-un 



for all A, B. For the resulting mass defect ratio one has 



Qvj-un 



exp ( - W (x lu ... uN )) = exp (-W ^ x^j j (3.31] 



and the simple linear superposition: 

N 



giu-ujy = %A- (3.32) 

A=l 



iii) For an aggregate of N fused particles of total bare mass M we obtain: 



A? 



l/i?iu...uiv = Yl Va/Ra =: , (3.33) 



A=l 

N 



fXA = ~w e (y0,1 ^ J2^ a = 1 - 

A=l 

iv) Using the asymptotic of function W(z) one easily obtains for the merger of N identical 
point particles, each of bare mass M and relativistic shift R 

( ( M\\ ln(iVf) 

giu ... UN = exp \-w [Njij) ~ -jjfHy when N ^°° ( 3 - 34 ) 

and 



R 1U ... UN = R, WN. (3.35) 

Here R is the relativistic shift of the separate particle. 

These formulas show that the accumulation of particles during the merger of N parti- 
cles leads to increase of the mass defect, since Qiu-un when iV — > oo - much like the 
situation in nuclear physics, as it should be from physical point of view. 



3.2.3 Some Solutions of the Mass Defect Equation for Continuous Distribu- 
tions 

To acquire some experience working with solutions of the MDE in the case of continuous 
mass distributions, we will consider in this subsection several simple examples of such 
distributions of identical particles with a constant density fi and a simple geometry: 

1. First we consider one-dimensional continuous mass distribution on a homogeneous 
circle of diameter d, and linear mass density n\, made of identical particles. 

From a symmetry reasons at the very circle one has (p = const. Then one has to 
calculate a simple integral in the rhs of Eq. (13.151) . Thus one obtains for the mass ratio 
the following expression: 

m - r ( w( - £ ki±£E]2*)) (3 . 36) 
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Figure 5: The dependence of the mass ratio Qdrde = m/M on the diameter d of circles 
with different fixed total bare masses M. The diameter d is shown in units of R. 



Suppose that the total bare mass of the circle is M. Then /ii = Mj nd. Replacing Hi 
in the relation (13.361) with this value, we obtain the result, shown in the Fig. [5j As seen, 
g(Q) = and when the masses of the circle are dispersed at bigger distances, increasing its 
diameter, the mass ratio increases and goes to 1 for d — > oo, as it should be from physical 
point of view. 

2. Our second example is a two-dimensional continuous mass distribution on a homo- 
geneous sphere of diameter d and surface mass density \ii = AM/nd 2 , made of identical 
particles. In this case one obtains in a similar way for the mass ratio: 

g{d) = exp ( - W (1 - j In (1 + d/R) )) Y (3.37) 

As seen in Fig. El when the masses of the sphere are dispersed at bigger distances, 
increasing the diameter, the mass ratio increases and goes to 1, when d — > oo, as it should 
be from physical point of view. 

In addition we see that in this example the mass ratio goes to zero for the non- 
physical value of sphere diameter d = —R. In contrast to the previous example, now 
one has g(0) = e - w ( 4M / R ) > o. This new property reflects the peculiar geometry of the 
spheres around the point sources: the limit of the sphere surface area remains finite when 
the sphere diameter goes to zero. 

3. For a homogeneous three-dimensional ball of radius r* one easily obtains the one 
dimensional integral equation: 

<p(r) = [ * dr'kiryy^ (3.38) 
Jo 

with kernel 



27T/1 . / . ., f r + r' + R * , 
k(r, r ) = r [r + r — |r — r | — Km I — ] J . (3.39) 
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Figure 6: The dependence of the mass ratio g sphere = m/M on the diameter d of spheres 
with different fixed total bare masses M. The diameter d is shown in units of R. 



Let us consider the limiting case of Eq. (I3.38j) with R — 0, and let us denote y?[ ] = v\r=o- 
Then we obtain the following Debye-Hiickel like boundary problem: 

^(0) = 0, £*^%)+^) = 0. (3.40) 

This problem has no auxiliary parameters and defines an universal function <P[o](£, £*). The 
solution gives the zero order term in the exact solution of the integral equation (13.381) : 
(p(r) = (p[o](y/ '4:7i fir) + O(R) and obviously describes the gravitational screening of the 
bare mass in the present simple case of three-dimensional continuous mass distribution. 
It is not hard to find the numerical solutions of the problem (13.401) . 



4 Some Concluding Remarks 

There exist at least three different approaches to the relation between physics and geom- 
etry: 

A) The classical physics considers the space-time continuum only as arena for the 
struggles of fields and particles. "These entities are foreign to geometry. They must be 
added to geometry to permit any physics" (see, for example, Misner and Wheeler in [IT]). 
In this approach the geometry is complete independent of physics. 

B) According to original Einstein's idea, the geometry of space-time is determined by 
the matter sources of gravity. Masses and non-gravitational fields are of non-geometrical 
origin. To some extend these sources of gravity are independent of geometry entities. In 
their presence the space-time continuum acquires the geometrical properties of curved 4D 
pseudo-Rieamannian manifold. Its geometry determines the motion and evolution of the 
very matter. Hence, in this approach we have a complicated interplay between geometry 
and physics, based on nonlinear differential equations. 



C) There is a third way, pioneered by Einstein and Rosen and developed in pure 
"geometrodynamics" by Wheeler and others. In this approach one may think that "There 
is nothing in the world except empty curved space. Matter, charge, electromagnetism 
and other fields are only manifestation of the bending of space. Physics is geometry", 
see Misner and Wheeler in |17j . In this approach we have "mass without mass, charge 
without charge", etc. At present the most well known hypothetical notions, created in 
the framework of this approach are the "black holes" and the "wormholes" of different 
type in space-time. They are consequences of the assumption to have " everywhere empty 
curved space" , see Wheeler in [T7] . 

In the present article we accept the original Einstein's idea, described in point B, about 
the relation between mater and geometry of space-time. Here we develop the mathemat- 
ical realization of this idea considering matter point particles as sources of gravity in GR. 
We have reached the following basic results: 

1. We have studied a new, two parameter class of solutions of Einstein equations. 
These static spherically symmetric solutions describe the gravitational field of massive 
point particle with bare mass M > and Keplerian mass m (0 < m < M). The difference 
between these masses, or their ratio g = m/M £ (0,1), defines the gravitational mass 
defect of the point particle. Such mass defect was not considered and studied until now, 
because for the standard Hilbert form (11.61) of the Schwarzschild solution "the bare rest- 
mass density is never even introduced" [27] . 

2. The new solutions form a two parameter family of metrics on singular manifolds 
M^ 1 ' 3 ^!^}, described in details in the present article, as well as in [T6l 122] . 

3. We have shown the principal role of the massive point source of gravity. Its presence 
offers a natural cutting for the physical values of the luminosity variable p £ [po,oo), 
where po > Pg- This happens because the infinite mass density of the matter point 
changes drastically the geometry of the space-time around it. This phenomenon is in 
sharp contrast to the situation in geometrodynamics, where luminosity variable may have 
an arbitrary close to zero value around the center of the black holes. 

A geometry of space-time with po = Pg > was discovered at first in the original 
article by Schwarzschild [3]. According to Eq. (I2.24p . such limiting value of the luminosity 
variable corresponds to zero value g = of mass defect ratio. For a finite value of m 
this is possible only if M = oo. In this sense our work is a proper extension of the 
Schwarzschild one to the physically and mathematically admissible values of the mass 
defect ratio g £ (0, 1). 

4. The existing attempts to describe the point matter source using Schwarzschild 
solution in Hilbert gauge (11.61) do not take into account the mass defect and thus fail to 
present a point idealization of the real relativistic objects. 

5. In full accord with Dirac's suggestion [2TJ our cutting of the domain of luminosity 
variable places the event horizon in the nonphysical domain of the variables. This effect is 
well known from the solutions of Einstein equations with massive matter sources of finite 
dimension [I]. 

6. The mathematical and the physical properties of the new solutions are essentially 
different in comparison with the well known other spherically symmetric static solutions 
to the Einstein equations. All of the new solutions have a strong singularity at the center 
of the symmetry, which is surrounded by empty space. They describe the single point 
particle sources of gravity in GR. The previously known static spherically symmetric 
solutions were often erroneously considered as a solutions, which describe a point mass, 



but this is not the case. 

7. It is clear that our solutions in generalized functions define in mathematical sense 
the fundamental solutions of the quasi-linear Einstein equations. These solutions are 
complete analogous to the fundamental solutions of Poisson equation in Newton theory 
of gravity. Thus the problem, formulated by Feynman in [TT] is solved. 

8. A proper quasi-linear superposition principle for initial conditions of Einstein equa- 
tions exist. It leads to a new theory of the relativistic gravitational defect of mass, 
illustrated in short in the present article. Our study shows that the hypotheses, used in 
this approach to the superposition principle in GR lead to physically reasonable conse- 
quences. One has to put on a more profound basis these hypotheses and to compare their 
consequences with the physical reality. 

9. Our results show that it may turn to be possible to transform the original Einstein 
relativistic theory of gravity, without change of its dynamical equations, into a normal 
physical theory with basic properties, which are intrinsic to the other branches of physics. 
This may permit us to get out of the way some of the specific scientific fictions, which are 
widespread at present and are thought to be an unavoidable consequences of GR. 

Especially, the correct theory of N matter particles systems seems to lead to the space- 
times with Euclidean topology, global time and standard physical causality. Hopefully, 
these important properties may transmute the now-days form of GR into a new version 
of the relativistic theory of gravity, compatible with the standard relativistic quantum 
mechanics of particles, much like the theories in the other branches of physics, and despite 
of the curvature of the space-time. 

10. Our basic conclusion is that we urgently need critical experiments and observa- 
tions, which can help us to choose the version of Einstein relativistic theory of gravity, 
compatible with the physical reality. Some idea of such type of experiments was recently 
proposed and discussed in [28J. 
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